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This is a mini-review of recent theoretical work in the field of relativistic heavy ion physics. 
The following topics are discussed: initial conditions and the Color Glass Condensate; ap- 
proach to thermalization and the hydrodynamical evolution; hard probes and the properties 
of the Quark-Gluon Plasma. Some of the unsolved problems and potentially promising direc- 
tions for future research are listed as well. 



1 Introduction 



In general, theorists get attracted to relativistic heavy ion physics because it is placed at the 
intersection of three different, and equally interesting, directions in contemporary theoretical 
research: i) small x, high parton density QCD; ii) non-equilibrium field theory; and iii) phase 
transitions in strongly interacting matter. Indeed, understanding the evolution of a heavy ion 
collision requires a working theory of initial conditions, of the subsequent evolution of the pro- 
duced partonic system, and of the phase transition(s) to the deconfined phase. This mini-review 
is an attempt to capture some of the recent changes and developments in the theoretical picture 
of these phenomena which have been triggered by an intense stream of the new data from RHIC. 



2 Initial conditions and global observables 

2. 1 The role of coherence 

Not so long ago, before the advent of RHIC, it was widely believed that at collider energies the 
total multiplicities will become dominated by hard incoherent processes. The very first data 
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Figure 1: Centrality dependence of the charged particle multiphcit y ne ar mid-rapidity in Au + Au collisions at 
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from RHIC (see^and references therein in this volume) provided a lot of food for new thought: 
the measured charged hadron multiplicities in Au — Au collisions appeared much smaller than 
expected on the basis of incoherent superposition of hard processes. Given that any inelastic 
rescatterings in the final state can only increase the multiplicity", we have an experimental proof 
of a high degree of coherence in multi-particle production in nuclear collisions at RHIC energies. 



2.2 Semi-classical QCD and hadron multiplicities 

Combining the idea of coherence with the parton model, we have to consider the initial parton 
wave functions of the colliding nuclei as coherent superpositions of the wave functions of the 
constituent nucleons. Since at small Bjorken x all of the partons in the nucleus at a fixed 
transverse coordinate participate in a hard scattering process, this treatment naturally leads to 
the notion of parton density in the transverse plane - a new dimensionful scale of the problem. 
Once this scale becomes comparable to the resolution scale determined by the kinematics of the 
hard scattering, the amplitude of the process is severely affected by the coherence. The limit 
on the parton density is reached when the occupation numbers of the gluon field modes with 
transverse momenta pT < Qs reach the value ~ l/as{Qs), characteristic for classical gauge 
fields - this is the phenomenon known as "parton saturation'El, leading to a coherent state of 
gluons - Color Glass Condensate (for reviews, see" 
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Since the integrated multiplicities are dominated by momenta pT < Qs and parton density 

1 /3 

in the transverse plane scales as ~ ^part (where Npart is the number of nucleons which 
participate in the process), Color Glass Condensate leads to a simple predictioiPfor the centrality 
dependence of hadron multiplicity in heavy ion collisions: 

duAA 



dr] 



Npart In (Npart) ■ (1) 



Combined with the dependence of the gluon structure function on Bjorken x known from HERA, 
which implies Q'^{x) ~ ^/x^, one can generalize this formula to predict the energy, centrality, 
rapidity, and atomic number dependencies of hadron multiplicitie:^ . Additional information on 
the dynamics of the collision can be inferred from the numerical lattice simulationi^ . So far 
this approach has been quite successful in predicting the multiplicities measured at RHIC; a 
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recent important example is given at Fig^ which shows the evolution of centrality dependence 
with energy in the entire RHIC range between ^/s = 20 and 200 GeV. One can see that the 
shape of the centrality dependence changes very little over a large energy range, in which the 
perturbative minijet cross section grows by over an order of magnitude. The prediction of the 
saturation model'"' is seen to agree with the data reasonably well; this indicates the possibility 
that parton saturation sets in in heavy ion collisions already at moderate energies. We do not 
expect the method to apply below y/s = 20 GeV however, since at lower energies the coherence 
length becomes shorter than the nuclear radius. 
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2.3 High px hadron suppression at forward rapidities, and quantum evolution in the Color 
Glass Condensate 

Parton saturation at transverse momenta < Qs sufficiently small x appears to have non- 
trivial consequences also for the nuclear dependence of the semi-hard processes. At very small x, 
when a.. In 1/x ~ 1, a semi-classical description has to be modified due to the quantum evolution. 
Small X evolution introduces anoma- 
lous dimension 7 ~ 1/2 in the gluon 
densities, so that the dependence 
on the momentum scale Q is mod- 
ified, to Q2 ^ g27^ 

Since in the 
vicinity of the saturation boundary 
the only dimensionful scale charac- 
terizing the system is the satura- 
tion momentum Ql, the cross sec- 
tion of semi-hard scattering should 
scale as a function of /Q"^ - it was 
foun d that this "geometrical scal- 
jj^g)i | ll | 12 | 13 | ll ] consistent with 

HERA data on deep-inelastic scat- 
tering. Combining these two ob- 
servations with the A dependence 
of the saturation momentum 

^1/3 

we come to the conclusio: 
that at sufficiently small x and mod- 
erate kT the nuclear dependence 
of hard processes in AA collisions 
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Figure 2: Nuclear modification factor in dAu collisions as a func- 
tion of transverse momentu m fo r different rapidities; the data from 
BRAHMS Collaboration, theoretical calculations fronff^. 



should change from SaQs " '■'•■part 
(where Sa ~ -^pirt the overlap 
area) to SaQ^P ~ ^-part- In 
(or dA) collisions the nuclear depen- 
dence is then SaQV ~ -4^/^ so 

there has to be a suppression as well. This suppression has also been found^''' in the numerica l 
solution of the Balitsky-Kovchegov equation, as well as i for recent work, see also' ' ' '. 

The experimental test of these ideas has been performed shortly afterwards - it has been 
established (for a review, see^^ in this volume) that at mid-rapidity y = there is no high 
kT suppression in dAu data; this means that the suppression observed in AuAu collisions has 
to come from the final-state effects, which will be discussed below. The data thus rule out the 
possibilit;P^ that x is small enough for quantum evolution to develop already at mid-rapidity 
at RHIC. Nevertheless, the presented arguments should apply at sufficiently small x. This is 
why the data on high kx hadron production at forward rapidities giving access to much smaller 
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Elliptic flow near mid-rapidity in Au + Au collisions as a fu nction of centrality (left) and transverse 

momentum (right); from ReiP^. 

values of x were eagerly awaited. The results from the BRAHMS experimenP^ demonstrated 
a strong suppression of high kx hadrons; moreover, the centrality dependence appeared consis- 

— 1/2 

tent with the predicted RdA ~ ^part scaling, in a dramatic contrast to the increasing Cronin 
enhancement observed at mid-rapiditv. Complementary results have been reported 
STAR Colaboration has also reportecP^ on an observation of a predicteJP^ nuclear-dependent 
weakening of the back-to-back correlations for hadrons separated by several units of rapidity. It 

will be interesting to check if the suppression extends to charm hadrons at forward rapiditiei?^!, 

on 

at mid-rapidity, the first results have been reported iix . Altern ativ e explanations based on 
"conventional" shadowing and multiple scattering (for a review, se^^ are also being explored. 



3 Approach to thermalization, and the role of classical fields 

There is by now an ample evidence of the importance of final state interactions in heavy ion 
collisions. Among the bulk observables, the azimuthal anisotropy of hadron production is a most 
spectacular evidence of this - indeed, if all of the elementary nucleon-nucleon collisions were 
independent, the produced hadrons would not be correlated with the nucleus-nucleus reaction 
plane . The observed azimuthal anisotropy (or the "elliptic flow", in the parlance of the field; see 
I 31 | 32 | fQj, review) indicates the existence of a correlation between the geometry of the nucleus- 
nucleus collision and the momenta of the emitted hadrons. An economical way of describing 
the evolution of a large number of particles in space and momentum is provided by relativistic 
hydrodynamics, which transforms the gradients of the initial parton density into the momentum 
flow of the produced hadrons. Hydrodynamical description is valid when the mean free path 
of partons is much smaller than the size of the system, i.e. when the system is sufficiently 
thermalized. The free expansion (or "inflation") of the produced system with time reduces the 
density gradients, so the magnitude of the elliptic flow crucially depends on the thermalization 
time when a hydrodynamical calc ulation is initiated. It appears that to describe the elliptic 
flow of the observed magnitudel^, one has to assu me that the thermalization time is very short, 
about Ttherm — 0.5 fm (for a review, see^^^^^^. Such a short thermalization time presents 
a problem both for the traditional perturbative and non-perturbative treatments. Indeed, in 
perturbative QCD the rescattering amplitudes are suppressed by powers of the coupling a^, 
so the thermalization time app ears long, on the order of 10 fm, which makes the description 
of the elliptic flow problematic In non-perturbative approaches, the interactions can be 
assumed strong, but the typical time scale of an interaction is ~ 1/Aqcd ~ 1 fm, so it is 
difficult to expect that several interactions needed for thermalization will occur during Ttherm — 



0.5 fm. The coherent classical fields present in the Color Glass Condensate (CGC) scenario 
may eventually provide a solution to this puzzle, since in this case the multi-gluon scattering 
amplitudes A{n — > m) from n ~ l/a^ to m ~ l/a^ glu ons are not suppressed. An approach to 
thermalization in this scenario was explored in Ref. r ecent 1 v. an attention was brought also 
to the role of instabilities in the equilibration process!"^. The use of CGC initial conditions of 
Ref. 1^ in a hydro dynamical approach!^ has led to a successful description of the RHIC data. 
Nevertheless, much work will have to be done to understand the thermalization process. 



4 Hydrodynamical evolution: more fluid than water 

Since hydrodynamical description relies on the direct use of the equation of state, the data can be 
used to extract an information on the properties of the medium. It app ears that the quark-gluon 
plasma equation of state as measured on the lattice (for a review, see"^ is successful in describing 
the data. However the data can tell even more about the prop ertie s of the medium, if one 
considers the influence of viscous corrections on various observables^^. Viscosity of the medium 
appe ars to affect the observables in a very significant way; in fact, one can deduce an upper 
limit^^^^on the ratio of shear viscosity r/ to the entropy density s, rj/s < 0.1 - much smaller 
than the same ratio for the water! Such a small value of viscosity, which reflects the dissipation 
of energy in a hydrodynamical evolution, contradicts the picture of weakly coupled quark-gluon 
plasma, and is more indicative of a strongly coupled quark-gluon liquid. A calculation of shear 
viscosity in the strong coupling regime of QCD i s st ill beyond the reach; however it has been 
made in = 4 super symmetric Yang-Mills theoryl^- the result is a small ratio of r//s = l/47r, 
comparable to the one inferred from RHIC data. 

A small value of viscosity in the strongly coupled quark-gluon plasma a posteriori justifies 
the use of the approach - to hadron multiplicities assuming the proportionality of the number 
of measured hadrons to the number of the initially produced partons. This assumption would 
be unnatural if the evolution of the plasma were accompanied by parton multiplication, but is 
justified if the viscosity is small and evolution of the system is close to isenthropic. 



5 High pt hadron suppression, jet quenching, and heavy quarks 

The suppression of high hadrons in Au — Au collisions is certainl y on e of the most spectacular 
new results at RHIC (for a comprehensive review of the data, seel^ in this volume) . Such an 
effect has not been seen at lower energie^; moreover, the results from the dAu run at RHIC 
indicate that at pseudo-rapidity ry = the observed suppression is entirely due to the final 
state effects, very likely a jet quenching in the quark-gluon plasma (for an overview, see'^^^^). 
Alternative scenarios, e.g. the absorption in a dense hadron gas, seem unlikely in view of the high 
energy density e ~ 20 GeV/fm^ (see e.g. achieved in the collisions. Nevertheless, additional 
experimental checks have to be performed; an important additional test of the jet quenching 
scenario involves the measurement of the suppression for heavy hadrons containing c or 6 quarks. 
If the suppression of high pt particles is indeed due to the induced radiation of gluons by fast 
partons, he avy quarks should lose significantly less energy than the light ones due to the "dead 
cone" effect'^. This prediction seems to be in accord with the first RHIC data, which within the 
error bar s in dicate no quenching effect on the spectra of open charm, as inferred from the decay 
electrons 1^^, however more precise data are desirabl e. Severa l other cal culatio ns of the energy 
loss of heavy partons have been performed (see Ref^^E^^^ and papers^^^^ in this volume); 
while they differ in the formalisms used, they all find a reduced energy loss for the heavy quarks. 
On the other hand, since heavy mesons {D,B, ...) have a typical large size determined by the 

"^A moderate amo unt of suppression in the SPS results however cannot be excluded due to uncertainties in the 
reference pp data^^ 



presence of the light quark in their wave functions, in the hadronic absorption mechanism one 
would expect that heavy mesons interact with about the same probability as the light ones. 
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in Au + Au collisions at = 200 GeV; from Re*^ 



6 Heavy quarkonium in hot QCD matter 

Ever since it was proposed as a signature of the quark-gluon plasmJ^, the dissociation of heavy 
quarkonia in hot QCD matter has remained a focal point of vigorous theoretical and experimental 
studies. The NA38/50 Collaborations at CERN have observed the suppression of J/ip and ijj' , 
and the current NA60 experiment will significantly extend the existing measurements (for an 



update on the recent res ults , se^^ in this volume). The first RHIC results have already been 
reported, in Au^JP^, dA'^^, and collisions. Theoretically, a new insight on the problem has 
been gained from the recent lattice calculations (sePEHIfor an overview) which indicate that the 
J/tp and rjc survive as bound states in the quark-gluon plasma at least up to the temperatures 
twice the critical, 2Tc. This observation is very important in understanding the properties of the 
strongly coupled quark-gluon plasma discussed above. However, in my opinion, it should not be 
interpreted as an indication that heavy quarkonia are not suppressed in the quark-gluon plasma 
unless the temperature is verv high; the point is that even if a cc state is bound in a plasma, 
it can be readily dissociated^- by the impact of gluons, which have much harder momentum 
distributions in a deconfined phasP^. Estimates of the activation rate of quarkonia due to the 
interaction with the heat batll^show that even if (cc) states exist as bound states, their yield 
can be strongly suppressed. More work has to be done to understand these effects better; on the 
lattice, a reliable extraction of the thermal widths of heavy quarkonia would be most desirable. 



7 Baryon dynamics 

The structure of baryons in non-perturbative QCD remains quite puzzling: while in non- 
relativistic quark model QQQ baryons are not so different from QQ mesons, in the approaches 
motivated by 1/Nc expansion they are drastically different - in Skyrmion picture, for example, 
they are the topological solitons of the meson fields. A closer look at the quark wave functions of 
baryons reveals that local gauge invariance requires the presence of novel configurations of gauge 
field - so called "baryon junctions"^"^. Naively, one expects that at high energies the collision 
of two relativistic nuclei would not lead to any substantial baryon stopping - since the valence 
quarks associated with the baryon number carry a large fraction of the nucleons' momentum, 
they are hard to stop in a soft process. However the account of non-perturbative baryon junc- 
tions leads to a substa ntial change in this picture, since the baryon number appears to be traced 
by soft gluon g"* in perturbation theory, baryon junctions were shown to correspond to 



multi-gluon exchanges in higher color representations' ■ . Substantial amount of baryon stopping, 

with the magnitude and rapidity dependence consistent with the baryon junction picture, has 

been observed at The influence of quantum evolution and parton saturation on the x 

distributions of valence quarks in nuclei has been addressed recently in Ref^^. 

Another exciting observation at RHIC related to baryon dvnamic s is a strong enhancement 

of baryon-to-pion ratios at moderate values of transever se m omentuml^ (so called -B/tt puzzle) 

and the larger magnitude of the elliptic flow for baryons^^. The proposed explanations include 

71 79 7*^ 

the phen ome non of parton coalescence^ - and the interplay of baryon junctions with jet 
quenchin^^. 



8 Summary 

The first years of the experiments at RHIC have changed in a dramatic way the theoretical 
picture of dense a nd h ot parton systems. The evidence for the existence of new states of QCD 
matter is mountin^^, and a consistent description of the observed phenomena has started to 
emerge. Nevertheless, hot and dense QCD is still in its infancy - and we have every reason to 
expect new surprises! 

I am grateful to the Organizers for the excellent meeting. This work was supported by the 
U.S. Department of Energy under Contract No. DE-AC02-98CH10886. 
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